Entanglement purification for high-dimensional multipartite systems 
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We propose an entanglement purification protocol for high-dimensional multipartite systems. In 
the protocol we can select a subensemble in a pure generalized Greenberger-Horne-Zeilinger (GHZ) 
state. This post-selection can be made by detecting the noise which contaminated the initial pure en- 
semble when the systems past through a noisy environment. For the detection of noise we investigate 
correlation properties of GHZ states and analyze their possible errors due to a noisy environment. 
We show that the presented protocol is more efficient than a simple generalization of the purification 
protocol for a bipartite state in high dimensions. 
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I. INTRODUCTION 

Entanglement lies at the heart of quantum informa- 
tion science. For instance, quantum entangled channels 
are essential in teleporting an unknown quantum state 
1] and in sharing a secret key for cryptography [2j. On 
the other hand, entanglement displays very fragile be- 
havior in a noisy environment: When a system initially 
in a pure entangled state comes into contact with an en- 
vironment, the state will decohere to a mixed state as 
it is unavoidably entangled with the environment and 
eventually will totally lose its entanglement [3[ . This de- 
coherence will decrease the degree of entanglement in a 
uantum channel for teleporation and degrade its fidelity 
Similarly, the decoherence will bring insecurity to 
the encrypted communication. Entanglement purifica- 
tion protocols (EPP's) have been developed to recover a 
subensemble in a maximally entangled pure state from a 
decohered ensemble before performin g re liable quantum 
information processes @, S 0, i, ©, EMl El 03 EI El . 
They have been applied to the distillation of Bell states of 
qubits d, [|| and many-qubit entangled states @, H, 1Ç | . 
An experimental realization has also been reported 11 1 . 



Recently, extensions of EPP's have been proposed for 
bipartite high dimensional systems and continuous vari- 
ables [H El El- 

Multipartite entangled states have been investigated 
intensively for quantum information processes, mostly 
based on the use of qubits. These studies in- 
clude Greenberger-Horne-Zeilinger (GHZ) nonlocality 
[ïï| , quantum secret sharing [Ï3| , cluster-state computa- 
tion Í18| . a nd quantum error correcting codes (QECC's) 
[Ï9l |20L l2l·l. |22| . Recently, strongly coupled systems have 
attracted interest due to their potential for scalable quan- 
tum computation and quantum simulation [23| . States 
of these systems are defined on high dimensional Hilbcrt 
spaces and are typically multipartite d-dimensional sys- 
tems, also known as qudits. It is desirable to understand 
what characteristics they have in terms of entanglement, 



what information processes they can be made to perform 
and if any advantages can be obtained over qubits. It 
has been shown that qudits enable more secure quantum 
cryptography than qubits 0, [H[. Two entangled qu- 
dits exhibit increased robustness against isotropic noises 
in nonlocality [26j]. Studies for many qudit systems in- 
clude nonlocality [27|, entanglement swapping Q, and 
preparation of GHZ states using hyper entanglement [11 . 
Optical realization of entanged qudits has recently been 
reported by using photon polarization I30l . orbital angu- 
lar momentum [3l|, and spatial modes [32|. It was also 
reported that a sequence of pump pulses leads to time-bin 
entanglement of qudits [11 . 

Entanglement purification is concerned with how 
to obtain a reduced number of copies of a high- 
fidelity entangled state from many polluted copies of 
noisy/nonmaximally entangled states, while quantum er- 
ror correction focuses on how to protect a certain quan- 
tum state from noise by encoding the state in a larger 
dimensional Hilbert space. Both approaches share a com- 
mon purpose of restoring the quantum coherence and 
they are closely related Q. Indeed, a QECC can be con- 
structed from a one-way EPP, where only uni-directional 
classical communication is allowed, and vice versa [yj. 
This fact was used for proving the security of quantum 
key distribution [H. EPP's may be classified into re- 
cursive, hash, and breed methods Recursive two- 
wayEPP's with bi-directional classical communication 
[EM S HL B EU can be constructed and analyzed in 
the frame of quantum error detecting codes (QEDC's) 
[3ol HE H3, while other EPP's utilizing local permuta- 
tions are also possible [II ■ 

In this paper, we propose a recursive two-way EPP for 
GHZ entangled states of many qudits. For the purpose 
we analyze correlation properties of generalized GHZ 
states (which will be called "GHZ states" throughout the 
paper). Based on the analysis, we find a measurement 
procedure to detect possible errors in the contaminated 
GHZ state as transmitted through a noisy environment. 
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Then, we present an EPP that can be employed for pu- 
rifing a many-qudit GHZ state. The presented protocol 
is shown to be more efficient than a simple generalization 
of modifying the protocol suggested previously to purify 
a two-qudit entangled state [l2| . 



by 



II. PERFECT CORRELATIONS FOR GHZ 
STATES 



Before considering the purification scheme, we first dis- 
cuss correlation properties of GHZ states. A maximal 
GHZ state shows perfect correlation properties for cer- 
tain composite observables and based on them it can be 
used for quantum communication. As subsystems in a 
maximal GHZ state are distributed through a noisy en- 
vironment, their interaction with the environment can 
contaminate the GHZ state with some errors and simul- 
taneously alters correlation properties in the GHZ state. 
By examining the alteration, we can reciprocally detect 
the errors caused by the noisy environment. Based on 
the detection result, we select copies with no errors. This 
post-selection increases the fidelity to the maximal GHZ 
state. Thus, investigating correlation properties of GHZ 
states are of importance in a purification scheme. 

Suppose that three observers, Alice, Bob, and Charlie 
want to share three qudits for quantum communication 
and further that the composite system of three qudits are 
in one of the GHZ states, defined by 
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where l,m,n € {0, 1, . . . , d — 1}, w = e l2 ^/ d , {\k)} is an 
orthonormal basis, and k Q m denotes k — m mod d, 
i.e., the residue when k — m is divided by modulus d. 
Here d is the dimension of Hilbert space Hd for a qudit 
and the subscripts A, B, and C stand for Alice, Bob, 
and Charlie, respectively. Note that the tensor product 
symbol "<£)" is omitted for a sake of simplicity in Eq. fï} . 
These GHZ states form a complete basis set in the Hilbert 
space Hd ® Hd ® Hd of the composite system. 

Each observer is allowed to choose one of two observ- 
ables, x and z. An observable is assumed to take, as 
its value, an element in the set S — {1, üj, . . . , uj^ 1 }. 
We introduce operators X and Z to represent orthogonal 
measurements for the observables x and z. An orthog- 
onal measurement is represented by a complete set of 
orthonormal eigenvectors and its outcomes are denoted 
by a set of eigenvalues. As the observables x and z are 
assumed to take an element in S, let S be the set of 
eigenvalues and then the observable operators are given 



X = J2u k \k) xx (k\, 

n=0 
d-1 



(2) 



n=0 



where \k) x ,z are eigenvectors of X and Z, respectively, 
corresponding to an eigenvalue u> k . The integer k of 
will be called an outcome throughout the paper. In this 
representation the observable operators X and Z are uni- 
tary 27] . To be more specific, let the measurement basis 
of Z be a Standard basis, that is, {\k) z = \k)}. The 
measurement basis of X is defined by quantum Fourier 
transformation F on {|fc)z}, i.e., each eigenvector of X 
is given by 



\k) x = F\k) í 



^ d-1 

Td^ 1 



(3) 



for an eigenvalue oj k . Then, the observable operator X 
is a level-shift operator in the Standard basis as 



d-1 



X = J2\ k ) ( k + 1 modd|. 



(4) 



k=0 



To describe a correlation, we employ a conditional 
probability among the outcomes p, q, and r that three 
observers obtain from their local measurements, respec- 
tively. Suppose Bob and Charlie try to infer Alice's out- 
come p, based on their obtained outcomes q and r. The 
conditional probability of their inference being correct is 
given by 



P(p\q,r) = 



P{p,q, r) 

P(q,r) : 



(5) 



where P(p, q, r) is a joint probability of Alice, Bob, and 
Charlie obtaining outcomes p, q, and r, respectively and 
P(q, r) is that of Bob and Charlie obtaining outcomes q 
and r. If P(p\q,r) ^ P (p), Alice's outcome is correlated 
to those of Bob and Charlie. The correlation is said to 
be perfect when a pair of observers, based on their out- 
comes, can infer with a certainty the other's outcome: 
P(p\q, r) = 1 for certain p, q, and r. 

To investigate correlations for a given GHZ state 
\ipimn) m Eq. ([T]), we need joint probabilities and con- 
ditional probabilities for composite measurements. A 
composite measurement is denoted as xaxbxc if three 
observers perform all X measurements on their received 
qudits, respectively. Similarly we denote another as zazb 
if Alice and Bob perform Z measurements with Char- 
lie performing no measurement. From the measurement 
xaxbxc-, the observers have joint probabilities given by 
Eqs. (P), (ÜJ), and quantum expectation, 



■ (p, Q, r) 



\xxx(p, q,r\ipi m n 1 1 

i 



d 2 



5{p- 



l mod d), (6) 
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where 5(k = mod d) is defined to be unity if k is con- 
gruent to zero modulo d and zero otherwise. Summing 
Px A x B x c ÍP> li r ) over Alice's outcome p leads to the joint 
probabilities P XgXC (q,r) of Bob and Charlic obtaining 
outcomes q and r: 

d-i ^ 
Px B xc(q,r) =^2Px A x B xc(P,q,r) = -^. (7) 

Note that the composite measurement xbxc, when Bob 
and Charlie both measure in X, with Alice performing no 
measurement, results in the same joint probabilities as 
Eq. |(7J), leaving no ambiguity in the expression ll P XBXC ." 

The conditional probabilities for the composite mea- 
surement xaXbXc are now obtained by substituting 
Eqs. © and Q into Eq. ©, 

Px A x B x c {p\li r ) = S(p + q + r = l mod d). (8) 

If p + q + r = l mod d, P XAXBXC (p\q,r) = 1, implying 
perfect correlations among the pertaining outcomes p, q, 
and r in xaxbxc- Following similar procedures for the 
given GHZ state \ipi mn )i we obtain conditional probabil- 
ities, 

p z A z B (p\q) = 5(p-q = m mod d), (9) 
Pz A z c (p\ r ) = S(p-r = n mod d), (10) 

for composite measurements za zb and z^zp, respec- 
tively. 

The perfect correlations for a GHZ state \il>imn) 
were shown in terms of the conditional probabilities in 
Eqs. (|8j)- (fTÜ|) . In an equivalent but alternative way, they 
can be characterized by showing that a set of secular 
equations satisfy with respect to some composite observ- 
able operators, called correlation operators. Consider 
three correlation operators, X ® X X 7 Z ® Z< ® 1, and 
Z ®1® Z* , with 1 the identity operator. These will be 
denoted as XaXbXc, ZaZ b , and ZaZ c , respectively, 
for a sake of simplicity. The perfect correlations revealed 
by Eqs. (f5))- (p~0]) correspond to the following cigcnvaluc 
relations: 

XAX B Xc\lpl,nn) = £*>'|Vw), (11) 
Z A Z B \i\)i mn ) =UJ m \lplmn), (12) 
Z A Z c \ip hnn ) =Ul n \lpl mn ). (13) 

It is clear that the subscripts Imn in GHZ states IV'Zmn) 
stand for the eigenvalues, or strictly their exponents, for 
the three correlation operators. The characteristics of 
perfect correlations in Eqs. ©-(0 or Eqs. |[ÏT |) -([Ï3"] ) play 
a crucial role in detecting transmission errors caused by 
the interaction with a noisy environment. 

III. DETECTION OF TRANSMISSION ERRORS 

When it is transmitted to a remote place, a system in- 
evitably interacts with an environment. The interaction 



causes undesired changes in the system. These changes 
may even include losing the system. To simplify our er- 
ror model we exclude such a possibility, which may be 
overcome by means of post-selection. The changes in 
a given state of the system are regarded as errors and 
are described by a trace-preserving completely positive 
quantum operation [38]. Let p be an initial state and 
p' be a contaminated state when transmitted through a 
noisy environment. The contaminated state p' can be 
expressed in Kraus representation by 

p' = Y;K a pKl (14) 

ct 

where Kraus operators satisfy the completeness relation, 
^2 a K\K a = 1. Each Kraus operator K a stands for an 
error during the transmission. We will introduce an op- 
erator basis, known as an error basis, to represent Kraus 
operators by linear combinations in the error basis. The 
error basis is assumed for a qudit to be 

E = {e i>j = Z i Xi\i,je{0,l,...,d-l}}, (15) 

where unitary operators X and Z are defined in Eq. © 
and are called level-shift and phase-shift errors respec- 
tively when used for representing errors. Note that the 
error basis E forms an orthogonal basis on the Hilbert- 
Schmidt space of operators as Tr(e| = d5i t kS^i 

where S^k is a Kronecker delta. 

The correlation properties discussed in Sec. |TT] play a 
crucial role in detecting transmission errors. To see this 
explicitly, consider the distribution of three qudits in a 
GHZ state IV'ooo) to three remote parties, Alice, Bob, 
and Charlie over a noisy environment. The initial state 
l^ooo) is an eigenstate with the same eigenvalues u> — 1 
for all the correlation operators, XaXbXc, ZaZ b , and 
ZaZç, as discussed in Sec.|Hl When errors oceur during 
the transmission, the initial state |^>ooo) changes. If er- 
rors belong to E, they change it to another GHZ state in 
Eq. |T]). It is obvious that the changed state can not have 
the same eigenvalues u>° for all the correlation operators. 
Such an alteration can be utilized in order to detect trans- 
mission errors. For example, supposing that a level-shift 
error X oceurs on Alice's qudit, then the state changes 
to Xa IV'ooo) = l^on)) which is an w-eigenstate of ZaZ b 
according to Eq. p^|) . If Alice and Bob measure Za 
and Zb on the contaminated state |^oii)) respectively, 
then their outcomes p and q satisfy p — q = 1 as implied 
by Eq. Due to the relation p = q being satisfied 

if no error oceurred, the error Xa can be detected in 
the composite measurement zazb- On the other hand, 
supposing that a phase-shift error Z oceurs on B, then 
the state changes to Zb \4>ooo) = l^ioo); which is an u>- 
eigenstate of XaXbXc as shown in Eq. (fTTj) . If Alice, 
Bob, and Charlie all measure X's on their respective qu- 
dits, then their outcomes p, q, and r satisfy p + q + r = 1 
according to Eq. ([5]) . The error Zb can thus be detected 
in the measurement xaxbxc as p + q + r = should 
be satisfied if there was no error. In a similar way, any 
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errors in E can be detected by appropriate composite 
measurements. Note that we have so far described a di- 
rect measurement on a single copy of a GHZ state for 
error detection. On the other hand, composite measure- 
ments must collectively be performed on several copies of 
GHZ states in QECC's/QEDC's [Ü^. This is also the 
case in our EPP, which will be presented in the following 
section. 

We have assumed that the errors occurring belong to 
the error basis E, neglecting general errors of linear com- 
binations in E. In this assumption, if an error in E oc- 
curs, the initial GHZ state |^>ooo) changes to another GHZ 
state. Due to stochastic occurrence of errors, the con- 
taminated state p' is a statistical mixture of GHZ states. 
In other words, p' is in a diagonal form in the basis of 
GHZ states, defined in Eq. fï}. O n the other hand, if 
more general errors are allowed, the contaminated state 
p' may also include off-diagonal elements in the GHZ 
basis. However, we use the measurements XaXqXc, 
ZaZb, and ZaZc for error detections in our EPP. As 
their basis are GHZ states (see Eqs. (fTTj) - ([TÏÏ]) ). their out- 
put states from p' will be statistical mixtures of GHZ 
states even if p' is not. Letting p" be a purified state 
from the input p', the diagonal elements of p" depend 
only on those of p' . If purified successfully, a certain di- 
agonal element of p" is close to unity. The independcncc 
of the off-diagonal elements in p' implies that neglecting 
general errors does not influencc the performance of our 
EPP. The concise model of errors that belong to E suf- 
fices for the error detections and significantly simplifies 
the analysis in our EPP. These kind of error models have 
also been employed in other EPP's @, 0|. 

IV. EPP FOR GHZ STATES 

Entanglement purification is concerned with obtaining 
a reduced number of copies of a high-fidelity entangled 
state from many polluted copies of noisy entangled states: 
It restores the quantum coherence of system polluted by 
environmental noise. In this section we propose an EPP 
for qudits which were initially in a GHZ state, say IV'ooo)- 
As they were transmitted in a noisy environment, the qu- 
dits have been affected by the stochastic errors modeled 
in Sec. IIHI and their state becomes a statistical mixture 
of GHZ states, given by 

P= Plmn \lplmn) {Í>lmn\ ■ (16) 

Z,m,n— 

The copies that have transmission errors can be detected 
by an appropriate composite measurement due to the 
perfect correlations of GHZ states, as discussed in Sec. [TT] 
and IIIII Those copies are discarded by means of post- 



selection based on the measurement result. Repeating 
the whole process, the state of copies remaining is close 
to the error-free state |V>ooo)) provided that the fidelity 
of p to IV'ooo) is higher than a certain threshold value. 
The fidelity pooo = (V'oool P |V>ooo) is used as a measure 
to investigate the performance of the EPP. If the fidelity 
Pooo converges to unity, then the purification is said to 
be successful. 

The present EPP is described in more detail for three 
parties sharing qudits in the mixed state p in Eq. (|16p . 
In this protocol every party performs scqucntially two 
subroutines PI and P2, that respectively correspond to 
the composite measurements xaxb%c and zaZbZc for 
detecting phase-shift and level-shift errors. They are 
schematically shown by two dashed boxes PI and P2 in 
Fig. [TJ PI and P2 are designed to eliminate phase-shift 
and level-shift errors in the error basis E (Eq. (JTSJ) ) and 
are performed sequentially. To eliminate phase-shift er- 
ror, the three parties receive two copies of GHZ states 
and all perform the subroutine PI on those. Alice, for 
instance, performs the inverse Fourier transformation 
on qudit 1 from one copy, the Fourier transformation F 
on qudit 2 from the other copy, and generalized XOR 
(GXOR) on them with qudits 1 and 2 as a control and 
target qudits, respectively. Here, GXOR is defined by 

GXOR\i) c \j) t = \i) c \iQj) t , (17) 

where the subscript c (í) indicates a control (target) qu- 
dit. Alice measures qudit 2 in the Standard basis (defined 
below Eq. ([TJ). The other parties do the same proce- 
dure on their qudits. After announcing their outcomes 
p, g, and r, the three parties keep the control qudits if 
p + q + r = mod d is satisfied. Otherwise, they discard 
those qudits. They repeat the previous process, until 
keeping two control copies. Applying the Fourier trans- 
formation F to the two copies completes the subroutine 
PI. The subroutine P2 starts by applying GXOR to the 
two copies kept, for instance Alice applies GXOR to qu- 
dits 1 and 3. Alice measures the target qudit 3 in the 
Standard basis, as in the box P2, Fig. [TJ The other par- 
ties do the same. All the parties announce their outcomes 
p, q, and r. If all the outcomes are the same, p = q = r, 
then they keep the control qudits, otherwise they discard 
them. The two post-selections will reduce the errors and 
increase the fidelity of the copies remaining. The whole 
process is repeated until a sufficiently high fidelity pooo 
is obtained. 

To analyze the performance of the present purifica- 
tion protocol, consider two copies in a composite state 
\ipimn) <S> l^/'m'n') which the parties share. In the sub- 
routine PI, just before the measurements, the composite 
state is transformed to 
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Bob 




FIG. 1: A schematic configuration of the entanglement pu- 
rification for three-qudit GHZ states. Alice, Bob, and Charlie 
share four noisy copies, labeled as 1,2,3 and 4. They perform 
local measurements and communicate via classical channcls 
to detect any transmission errors. Dashed box PI (P2) is de- 
signed to detect and eliminate phase-shift (level-shift) errors. 
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p'+g'+r' = — l· 



Ojímq+nr-m'q'-n'r') ^ ^ r )\p Q p' , q Q q' , r Q r') , 



(18) 



r 



where "p = V implies "p = l mod d." After their mea- 
surements on the target qudits, the parties sum the out- 
comes and obtain {p Q p') + (q 8 q') + (r 6 r') . The sum 
is equal to l + V mod d, noting that p + q + r = l mod d 



and p'- 



= —V mod d. If the sum l+l' = mod d, 



the parties keep and perform F on their control qudits, 
which are now in the state \ipi,mQm',nen')- This complete 
the subroutine PI. In the subroutine P2, consider two 
copies in a composite state \ipimn) 8> iV'í'm'n') which are 
kept in PI. In their measurements on the target qudits 
following the GXOR operations, the parties will have all 
the same outcomes if m — m' and n = n' , and keep the 
control qudits. In this case, these qudits are in the state 
\^>i®i' ,m,n) where l © V denotes l + V mod d. 

When two input states in PI are statistical mixtures 
P = J^Pirnn \ipimn) (ipimn \ in Eq. ([Tïï]) , the output state is 
a new statistical mixture p' = J2i. m , n Pi m n \^lmn) (4>lmn\ 
with recurrence relations between pi mn and p\ mn ■ One of 
these is given by 



Pooo 



mnPOmn 



El 



,771,71,771 .71 



-Q Plvn7iPdQl,m' n f 



(19) 



The output state in P2 for the given input p 1 becomes 
another mixture p" = Y,i, m ,nPlmn\^lmn) {^lmn\- The 
recurrence relation between p'q QQ and p\ mn is given by 



Pooo 



Ed— 1 / / 
1=0 PiooPd 



ei, oo 



1,1' ,m,n=0 Plr, 



(20) 



iPl' 



If póoo > -Pooo an d Pooo > PÓoo; then we obtain a purified 
state through the subroutines PI and P2. The condi- 
tions are reduced to a threshold fidelity such that if an 
initial fidelity is higher than it, the purification will suc- 
ceed. With the help of the recurrence relations includ- 
ing Eqs. (|19|) and (|20p . we can numerically determine 
whether the present protocol successfully extracts |"0ooo) 
from statistical mixtures. The present protocol will be 
called a direct EPP. 

We shall now present another protocol by generalizing 
the proposals of Ref. [Ï2j for high-dimensional Bell states 
of qudits. This protocol is an indirect EPP in the sense 
that GHZ states will be purified through the purification 
of Bell states in the way. It may be regarded as a simple 
extension from the bipartite EPP, as discussed by other 
authors [7]. The indirect protocol takes the following 
steps: 

(1) An ensemble of qudits in noisy GHZ states, p 
in Eq. (|16p . is divided into an equal amount of two 
subensembles. 

(2) From each subensemble, Alice, Bob and Charlie 
share a trio of qudits in a noisy GHZ state. For the 
trio of qudits from one subensemble, Bob measures his 
qudit in the basis {1^)^} and announces the outeome k 
to the others. Receiving the outeome k, Alice applies an 
operation Z k to her qudit, where Z k is the k-th power of 
Z in Eq. ((2|). Then, Alice and Charlie have the pair of 
qudits in a reduced (high-dimensional) Bell state. From 
the other subensemble, Charlie and Alice do the same as 
Bob and Alice did, and Alice and Bob have another pair 
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TABLE I: Threshold fidelities /dir and ,/md, the minimum 
fidelities required for the successful purification by the direct 
and indirect EPP's are presented for noisy GHZ states, as 
given in Eq. (|2ip . The number d is a dimension of the Hilbert 
space for each qudit. 

d /dir /ind 



2 0.4073 0.4167 

3 0.2305 0.3000 

4 0.1555 0.2227 

5 0.1155 0.1780 

6 0.0907 0.1489 



in a reduced Bell state. This step is repeated over the 
subensembles. Alice and Bob have one subensemble of 
qudits in noisy Bell states p\ and Alice and Charlie have 
the other in pi. 

(3) Alice and Bob (Alice and Charlie) purify their noisy 
Bell states, p\ (/O2), by performing the EPP by Alber 
et. al. 12]. They obtain a reduced number of copies in 
purifred Bell states, p[ (p'2), if p has a fidelity higher than 
a certain threshold. 

(4) From each post-selected subensemble, Alice 
chooses a pair of qudits. She then share one with Bob and 
the other with Charlie. She performs a GXOR operation 
on her qudits, where a control (target) qudit belongs to 
the pair shared with Bob (Charlie). Measuring the tar- 
get qudit in the Standard basis {|&)}, she announces the 
outcome k and Charlie applies X dok to his qudit. They 
obtain then a state closer to |"0ooo) ■ This step is repeated 
over the post-selected ensembles. 

We compare the direct and indirect EPP's in terms 
of their threshold fidelities, fdh and /i n( j respectively, 
for successful purification and the minimum numbers of 
copies required to achieve a certain high fidelity, say 0.99. 
For the comparison, consider a mixed state of 

a#ooo)<^ooo| + -^3^1, (21) 

where t/d 3 is a purely random state of three qudits. In 
Table U we present the threshold fidelities /di r and /i n d 
in terms of dimension d. The thresholds decrease as the 
dimension d increases, implying that both EPP's work 
for lower fidelities of mixtures at higher dimensions. For 
each dimension d, the direct protocol has a lower thresh- 
old than its counterpart and in a way the direct EPP 
can be seen to work better than the indirect. The direct 
protocol is also more efficient in the sense that it requires 
a smaller number of copies to achieve a certain high fi- 
delity. For instance, considering a case of d = 6, suppose 
that the initial state in Eq. (|2"Tj) is of fidelity 0.5. In the 
direct protocol, about 48 copies suffice to achieve fidelity 
0.99. On the other hand, the indirect EPP needs about 
192 copies. These comparisons in terms of the threshold 



fidelities and the minimum numbers of copies leads us 
to conclude that the direct protocol shows better perfor- 
mance than its counterpart. It is because the direct pro- 
tocol utilizes multipartite correlations, whereas the in- 
direct one employs only bipartite correlations. Utilizing 
multipartite correlations leads to efficient entanglement 
purification. 



V. REMARKS 

With the present direct EPP we perform alternatively 
the subroutines PI and P2. Other combinations of sub- 
routines, for instance P1+P1+P2 are potentially possi- 
ble. Such variants may be beneficial against some trans- 
mission errors. For instance, it was argued in Ref. [39[ 
that phase errors are dominant in experiments for high- 
dimensional entangled states. In this case the subroutine 
PI could be employed more frequently than P2. A de- 
tailed analysis of such variants is beyond the scope of this 
paper. 

A purification protocol utilizes a detecting scheme of 
errors rather than a correcting scheme. However, the 
results from QECC's are also important as they can 
easily be adopted to QEDC's. QECC's, in particular 
Calderbank-Shor-Steane codes treat separately bit-flip 
and phase-flip errors [ï{|. This treatment of separat- 
ing errors was shown to be useful in developing EPP's 
for qubits 0], which have in common two subroutines 
similar to our PI and P2. 

In summary, we have investigated and analyzed the 
correlation properties of GHZ states for qudits (d- 
dimensional systems) . It was shown that a noisy environ- 
ment alters the correlations of a GHZ state if it contami- 
nates the state of qudits during the transmission. Bascd 
on this, we have suggested a detection scheme for errors 
and then a direct EPP that purifies noisy GHZ states 
of qudits. By comparing with the indirect EPP derived 
from the EPP's for two-qudit Bell states, we showed that 
the direct EPP shows better performance than the indi- 
rect one in terms of threshold fidelities and the minimum 
numbers of copies for achieving a high fidelity. 
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